Numerical Calculation of the Neutral Fermion Gap at v — 5/2 
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The p = 5/2 fractional quantum Hall state Oj-yl] has been 
the subject of intense experimental and theoretical investi- 
gation in recent years because it may support non-Abelian 
anyons and may serve as a platform for topological quantum 
information processing 14j-|6[]. Theoretical 171 4 1311 and exper- 
imental 1 3-3] evidence has been rapidly accumulating in 
favor of v = 5/2 being an Ising-type non-Abelian state, in 
the universality class of either the Moore-Read (MR) Pfaffian 
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We present the first numerical computation of the neutral fermion gap, A,/,, in the v — 5/2 quantum Hall 
state, which is analogous to the energy gap for a Bogoliubov-de Gennes quasiparticle in a superconductor. We 
find ~ 0.027^-, comparable to the charge gap, and discuss the implications for topological quantum 
information processing. We also deduce an effective Fermi velocity vf for neutral fermions from the low- 
energy spectra for odd numbers of electrons, and thereby obtain a correlation length ^ = vf/A^ « 1.3 £o. 
We comment on the implications of our results for electronic mechanisms of superconductivity more generally. 
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fermions ip correspond to Bogoliubov-de Gennes quasiparti- 
cles in the superconductor. In most superconductors, these 
two gaps have completely different scales and are not consid- 
ered on the same footing. However, in the v = 5/2 state, there 
is only a single energy scale e 2 /e£o, so these gaps can be com- 
parable. Thus far, however, only A a has been computed. In 
this paper, we compute A^,. This is the appropriate quantity 
to use when comparing the gap in the v = 5/2 state to the 
gaps in other superconductors, and when drawing lessons for 
non-phonon mechanisms of superconductivity from this state. 

The neutral fermion gap is also a relevant quantity in de- 
termining the effectiveness of topological protection in the 
v = 5/2 state. The transfer of Ising ip charge between quasi- 
particles, e.g. through tunneling, alters the non-local state 
shared by the two quasiparticles. It is, thus, responsible for 
splitting the degenerate non-local states and causing errors in 
the encoded information [24]. Similarly, the neutral fermion 
gap directly determines the visibility of non-Abelian statisti- 
cal signatures in interference experiments (see Il25ll . and ref- 
erences therein), since tunneling of the neutral ip charge (be- 
tween bulk quasiparticles and between bulk quasiparticles and 
the edge) suppresses interference terms. In this light, it is of 
paramount importance to study this quantity. In this letter, we 
produce numerical estimates of the neutral fermion gap and 
correlation length for the v = 5/2 non-Abelian quantum Hall 
state. 

In order to model the v = 5/2 state, we assume that both 
spins of the lowest Landau level are filled and inert and focus 
on the second Landau level, which has v = 1/2. Our calcula- 
tion neglects finite layer-thickness 11311 and Landau-level mix- 
ing lid 11 , 2r3l . which certainly play a role in real devices. A 



18[] or the anti-Pfaffian (Pf) state GHHH- 
The potential use of this state for topological quantum in- 
formation processing is dependent on the size of the energy 
gaps A a to different species of quasiparticles a. If the temper- 
ature T can be kept much less than these gaps A a and inter- 
quasiparticle distances x kept much greater than the tunneling 
correlation lengths £ a , then the corresponding error rates will 
vanish as e _Aa / T and e~ x ^ a and, hence, be negligible. 

The smallest gap for charged quasiparticles is usually as- 
sumed to correspond to the minimally charged excitations of 
a state [37]. For the MR and Pf states, the minimal charge 
±e/4 quasiparticles also carry non-Abelian Ising topological 
charge a. It is natural to interpret the gap corresponding to the 
temperature dependence of the longitudinal resistance, p xx ~ 
e -A tran5 /2T^ as tne ener gy g a p ^qh _|_ ^qp f or a charge ±e/4 

quasihole-quasiparticle pair, which is thereby deduced from 
experiments to be A^ h +A^ p = A trans w 0.5K in the highest- 
mobility samples 12111 . Numerical studies of small numbers of 
electrons interacting through Coulomb interactions in the sec- 
ond Landau level find A^ h + Ag. p w 0.025 - 0.029^ (which 
is 3.2-3.7Kat6.5T) J10]- 



However, bulk electrical transport is not sensitive to the en- 
ergy gap of electrically neutral excitations, such as the neutral 
fermion that carries Ising topological charge ip in the MR and 
Pf states. Consequently, A^ has not been measured (though 
it could, in principle, be determined from thermal transport 
measurements or, as we discuss below, from interferometry 
measurements in mesoscopic devices). A^ has previously not 
been theoretically calculated, either. 

The MR and Pf states are the quantum Hall analogues of 
spin-polarized p x + ip y superconductors 03 01 El- Charge 
e/4 quasiparticles a correspond to flux hc/2e vortices; neutral 



more realistic calculation, including these effects, will be dis- 
cussed elsewhere. Here, we focus on the simplified situation 
of an infinitely-thin two-dimensional layer in a very high mag- 
netic field and study small systems (N e < 15 electrons) by 
exact diagonalization and larger systems (13 < N e < 26elec- 
trons) b y th e density-matrix renormalization group (DMRG), 
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On the sphere, the MR ground-state occurs at N^, = 2N e — 
3, where the number of electrons N„ should be taken to be 
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even. For the following discussion, it will be helpful to intro- 
duce the following terminology. For a given, fixed arbitrary 
electron number N e and flux N^, we will call the state of 
lowest energy the "lowest energy state." If N e is even and 
Nef, = 2N e — 3 we will call the lowest energy state the MR 
ground state. The reason for this distinction is that for some 
values of N e , N^, the lowest energy state should be under- 
stood as a state with a quasiparticle excitation. We denote the 
lowest energy for a system of N e electrons and fluxes by 
E(N^N e ). 

In order to compute the energy gap for an electrically- 
neutral quasiparticle, we need to compare the usual ground- 
state to a configuration that forces the system to have a neu- 
tral excitation of non-trivial topological charge. In an Ising- 
type system, the lowest energy state on the sphere with N e 
odd electrons must have non-trivial quasiparticles whose total 
topological charge is ip ll36ll . The two simplest possibilities 
are that such a state either has a neutral tp quasiparticle, or 
a charge e/4 a quasihole and — e/4 a quasiparticle pair that 
fuses into a ip. (If such a pair forms a bound state, it is equiva- 
lent to a neutral ip quasiparticle.) Which of these possibilities 
actually occurs depends on whether the energy to create 
a neutral fermion ip is less than or greater than the energy 
A^ 1 + Af to create the quasihole-quasiparticle pair. 

Consequently, E (N<j, + 2, N e + 1) - E (N4, N e ) is the 
energy £ due to one electron plus/minus the energy of the 
non-trivial quasiparticle(s) for N e even/odd (either or 
A| h + Act P ). To isolate the energy of these (collectively) neu- 
tral quasiparticle(s) with total topological charge ip, we define 
the neutral fermion gap 



A F (N e ) 
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[E (Ntj, + 2,N e + 1) 



-E(N <t> -2,N e -l)-2E(N <j> ,N e )] (1) 



for paired states. In the regime in which E(N^,,N e ) 
scales linearly with N e , the neutral fermion gap Ap (N e ) 
will be constant. It is instructive to contrast Eq. Q] 
with the expression for the charge gap, A c (N e ) = 
\ [E (JV + 1, N e ) + E{N^- 1, JV C ) - 2E (N^, N e )}. In 
Eq. Q] we compare the energies of systems with the same 
charge-flux relation so that the net charge of all excitations 
is zero while A c compares the energies of states with fluxes 
offset by one so that the net charge of all excitations is ve. 

For pure Coulomb interactions in the second Landau level, 
we have computed the ground state energies for even numbers 
of electrons up to N e = 26 and the lowest state energies for 
odd numbers of electrons up to N e = 17. In a recent cal- 
culation, Lu et al. 12811 have computed these energies up to 
N e = 18 electrons by exact diagonalization; our energies are 
in agreement with theirs. In Fig. Q] we show the values of 
the neutral fermion gap Ap(N e ), computed using Eq. Q] as 
a function of inverse system size 1/N e for up to N e = 17 
electrons. As may be seen from Fig. [U the neutral fermion 
gap fluctuates considerably, which is a sign of finite-size ef- 
fects. If we were to use a purely linear fit, then we would find 




FIG. 1: The neutral fermion gap, as defined in Eq. (TJ, for the v 
5/2 MR state as a function of inverse system size 1/N e . 



limjv e _>-oo Ap(N e ) w 0.028; if we were to fit the gap to a 
constant, we would find limiv e ->oo Ap(N e ) sa 0.023. How- 
ever, the errors in these fits, determined from the maximum 
fluctuation away from the average, are large (though Af is 
clearly non-zero). Therefore, more care is needed in order to 
perform an N e — > 00 extrapolation. 

To this end, we note that if the system is gapped, then we 
can write E (N^, N e ) in the form 



E{N^N e ) = £N e + E t 



even, odd 



0{e 



N, 



) (2) 



for N e even or odd, respectively. The leading terms are the 
same for even and odd N e because the energy per particle £ 
must be the same in the thermodynamic limit. The constant 
terms -E e ven,odd are due to the internal order of the phase and 
the genus of the system 03611 . as well as the energy cost of the 
(collectively) neutral quasiparticle(s) for N e odd. Corrections 
to these first two terms are exponentially small in the linear 
size of the system (~ y/Ne) since the system has a gap; here, 
a is a constant inversely proportional to the correlation length. 
Substituting Eq.|2]into Eq.Q] we find 



Af (N e ) = £U - E,. 



0{e 



(3) 



and thus limjv e -).oo AF(N e ) = E dd — E even , further justify- 
ing our definition of the neutral fermion gap. We can, how- 
ever, use Eq.|2]to extract E^dd ~ -Seven more directly by simply 
fitting the numerical data with functions of this form, and it 
allows us to exploit the larger system sizes for which we have 
computed the ground state energies for even N e . In Fig. [2] 
we plot E N e ) /N e vs. 1/N e , and fitting to Eq. |2] (di- 
vided by N e ) but replacing, for simplicity, the 0(e~ a ^^) 
term by a single term ce _av ^, we find £ = —0.3634, 
-Seven = -0.5381, and E oM = -0.5114. For N e even, we 
find c = —0.7876 and a = 0.6675, while for N e odd we find 
c = -1.4700 and a = 0.8287. Thus, we can reliably extract 
the thermodynamic limit of the neutral fermion gap by taking 
the difference between the 1/N e terms in the expressions for 
E (Nfi, N e ) /N e . We find E odd - E even w 0.027 (in units of 
e 2 /e£ ). 

One advantage of using this method of extracting the neu- 
tral fermion gap is that it is easier to diagnose potential dif- 
ficulties with the N e — > oo extrapolation. For instance, 
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FIG. 2: The N e — ► oo extrapolation of E {N^, N e ) /N e correspond- 
ing to the v = 5/2 MR state for iV e even (squares) and N e odd (dots) 
using Eq.[2] from which we find A,/, ~ 0.027. 



one potential pitfall is aliasing. If one of the systems stud- 
ied is actually in the ground state of a different phase, then 
E (Nj,, N e ) /N e would not sit on the expected (nearly-linear) 
curve. As may be seen from the figure, the data points deviate 
negligibly from the fitting curves, so this is not the case for the 
system sizes we study. At any rate, the most serious potential 
aliases occur at N e < 10, which we do not consider for the 
extrapolation in Fig. [2] 

Having computed lini7v c ^co Aj?(7V e ) 0.027 for the 
(thermodynamic limit of the) neutral fermion gap, we now 
address the nature of the associated quasiparticle. The state- 
ment that the v = 5/2 state is an Ising-type topological state 
merely guarantees that ip is an allowed value for the topologi- 
cal charge in any region bounded by a closed curve. It does not 
guarantee that there is actually an energetically stable quasi- 
particle that has this value of topological charge. The mis- 
match between allowed topological charges and stable quasi- 
particle species is is a feature of all to polo gical states. For 
instance, in the v = 1/3 Laughlin state [29], the charge 2e/3 
quasihole carries an allowed value of topological charge, but 
it is not an energetically stable excitation (for Coulomb inter- 
actions); if we attempt to create one, it will decay into two 
charge e/3 quasiholes. Similarly, we must consider the pos- 
sibility that a neutral %l> quasiparticle will simply decay into 
a charge ±e/4 a quasihole-quasiparticle pair that fuses into 
the ip channel, i.e. that A^ > A^ p + Af 1 . In this case, 
lirriAr_._j.oo Ap(N e ) = _E 0( jd — E ewen would be identified with 
A^ p + A^ h and provides a lower bound for A^. However, 
since we find E dd — E even w 0.027 and previous studies II 1211 
obtained A^ p + A^ h w 0.029, we tentatively conclude that 
the neutral fermion tp is stable and has 



A, 



0.027. 



(4) 



Stronger evidence supporting this interpretation comes from 
the good fit of our data to the N e odd case of Eq. [2] If the 
ip quasiparticle were unstable, there would be a — 
term in the odd electron number energy, resulting from the 
Coulomb interaction energy between the ±e/4 charges ll_2Tl . 

For purposes of comparison, we note that a similar com- 
putation of the neutral fermion gap for the v = 1/3 Laugh- 



lin state 12911 would give the value zero because the even- 
and odd-electron number ground state energies lie on the 
same line 11211 : since it is not a paired state, there is no 
qualitative difference between even and odd electron num- 
bers. On the other hand, the Pf state, the (3J3, 1) state Q 
and the Bonderson-Slingerland (BS) states Bill have neutral 
fermionic excitations whose gaps can be computed by the 
method explained in this paper. In the absence of Landau- 
level mixing, A,/, is expected to be precisely the same for the 
Pf state as it is for the MR state; preliminary calculations are 
consistent with this, as we report elsewhere 13211 . In the case 
of the k > 3 Read-Rezayi states plsll . there are neutral excita- 
tions that are non-Abelian and, therefore, cannot be obtained 
by simply altering the electron number and flux. 

Although the neutral fermion gap has not been previously 
calculated, a related quantity has recently been calculated, 
namely the splitting between the two degenerate states that 
occur for four e/4 a quasiparticles 13411 . This splitting, 
AE(r), decays with distance r between the a quasiparticles 
as AE(r) ~ f(r) e~ r /'» for large r. Here, f(r) is an oscil- 
latory function and £ is the characteristic length scale for the 
decay. If we interpret this splitting as the energy associated 
with inter-quasiparticle tunneling of neutral fermions, then we 
expect £ = £/, = v/A^, where v is the velocity of a neutral 
fermion. If the MR state is interpreted as a paired state with 
small gap, then v would be the Fermi velocity vp of the un- 
derlying Fermi-liquid-like metallic state. In such a case, the 
Fermi velocity could be deduced by studying the spectrum of 
a single neutral fermion as follows. For odd N e , the energy 
spectrum will not have a gap above the lowest energy state 
(in the thermodynamic limit) since there will be one unpaired 
neutral fermion above the Fermi energy, and this fermion can 
be excited to any other state above the Fermi energy. In a 
BCS mean-field theory, the energy spectrum for odd N e will 

be bounded below by the curve El = 



E, 



2A ; 



■ Ap + E gs , where E gs is the ground state energy for 



N e — 1 electrons, is a single-particle energy relative to the 
Fermi energy for a state with angular momentum L. We take 
(l = j^[L(L + 1) - L (L + 1)), where L is the highest 
occupied angular momentum orbital. Thus, for L w L , the 
excitation energies are expected to be quadratic in L — Lq: 



E l 



1 ( v F (2L + 
2A F V Nfto 



1) 



(L - L a f + const. (5) 



As may be seen in Fig. [3] the lowest excitation energies for 
N e = 9,11,13,15 appear to follow a parabola. A linear 
extrapolation of the vf values obtained from these spectra 
according to Eq. [5] gives vf ~ 0.021 e 2 /e, which leads to 
£0 w 0.8 £q. However, the parabolic fit is quite poor for 
N = 13; the other three system sizes are consistent with 
vf ~ 0.035 e 2 /e, or £^ w 1.3^o- We note, for compari- 
son, that Baraban et al. I34I1 find a length scale £ sa 2.3^o, 
although their calculation is for much larger system sizes and 
for trial wavefunctions, rather than the Coulomb ground state. 
Our results imply that a quantum computer based on the 
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FIG. 3: The low-energy spectrum of the v = 5/2 MR state for 
9, 11, 13, 15 electrons with parabolic fits of the lowest-lying states 
to Eq.[5] from which we extract the effective velocities plotted in the 
lower panel. 



v = 5/2 fractional quantum Hall state should be operated at 
temperatures much lower than Ay,, which is rj 3.4K for a 
magnetic field B = 6.5T. This implies that, at 35mK, the rate 
at which phase errors can be expected for a topological qubit 
is ~ e ~ A ^/ T py 10~ 44 if the computational anyons are fur- 
ther than £y w 130A from each other and the edge. In the 
experiments of Willett et al. HH], the inter-quasiparticle dis- 
tances are probably comparable to this implies that the er- 
ror rate may be large and there is probably significant splitting 
between the 2™ _1 states expected for 2n quasiparticles 13511 . 
By measuring the time over which the signal through an inter- 
ferometer remains stable, it should be possible to measure the 
error rate and, thereby, Ay, . In addition, bulk thermal transport 
may be dominated by thermally-excited neutral fermions. Al- 
though charge e/4 quasiparticles may have a smaller energy 
gap (approximately half that for a ip), they will be much more 
strongly localized by disorder than neutral fermions. 

2 

Finally, we note that Ay, sa 0.027 ^- is small compared 



to the Coulomb energy. Since one might argue that the gap 
is small because of the proximity to competing phases, such 
as the striped phase |8[], we consider the neutral fermion gap 
for a Hamiltonian in which the only interaction is the (re- 
pulsive) three-body interaction for which the MR wavefunc- 
tions are the exact ground states H, 18]. For this Hamilto- 
nian, the ground state energy is precisely zero for N e even, 
so -E e ven = 0. Thus, we must only compute the ground state 
energies for iV e odd. A simple linear extrapolation of these 
energies gives lim.jv e ^cx3 Aj? = Edd « 0.45, if the coeffi- 
cient of the three -body interaction is 1. Thus, there is nothing 
wrong in principle with the naive idea that the superconduct- 



ing gap can be comparable to the Coulomb energy scale for an 
electronic pairing mechanism, so long as there are no nearby 
competing phases to suppress it. 
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